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1 Introduction 



The method of QCD sum rule |]T] has become a powerful tool in studying hadron physics on 
the basis of QCD. It is a framework which connects physical parameters of hadrons to the 
QCD parameters. In this approach hadrons are represented by their interpolating quark 
currents taken at large virtualities and then correlation function of these quark currents is 
introduced. Following this, on one side, the correlation function is calculated by the operator 
product expansion (OPE), and on the other side, its phenomenological part is constructed. 
Physical quantities of interest are determined by matching these two descriptions at large 
virtualities via dispersion relations. 

The QCD sum rule method is discussed in many review articles (see [@]-|| and references 
therein) and successfully applied in studying various characteristics of hadrons. One of the 
important characteristic static parameters of hadrons is their magnetic moments. Magnetic 
moments of nucleon were calculated in framework of the QCD sum rule approach in |7|, |§ 
using external field technique. They were later refined and extended to the entire baryon 



octet U 10 



In this work, we present an independent calculation of the magnetic moments of the 
octet baryons in the framework of an alternative approach to the traditional sum rules 
method, i.e., the light cone QCD sum rules method (LCQSR). In the LCQSR method OPE 
is carried out near the light cone x 2 = instead of at the short distance x ~ 0. This is an 
expansion over the twists of the operators rather than the dimensions, as is the case in the 
traditional QCD sum rules. The nonperturbative dynamics is parametrized by the so called 
light cone wave functions, rather than the vacuum condensates used in the traditional QCD 



sum rules. A detailed description of this method can be found in || [II], |12[ . There are 
many applications of the LCQSR in the current literature. Note that magnetic moments 
of the nucleon and decuplet baryons were studied in |H| and JR], respectively, using the 
LCQSR method. 

The paper is organized as follows. In section 2 the LCQSR for the magnetic moments 
of the octet baryons are derived. In section 3 we present our numerical calculations and a 
comparison of the prediction of the LCQSR method on the magnetic moments of the octet 
baryons with the results of the other methods and experimental results. 



2 LCQSR for the magnetic moments of the octet baryons 

As we have already noted, in order to construct the sum rules for the magnetic moments of 
hadrons, it is necessary to introduce the correlation function of interpolating quark currents 
which have the same quantum numbers as corresponding hadrons. For this aim we consider 
the following correlator function: 

n = ij d 4 xe ipx (0|T{77B(z)rjB(0)}|0) 7 , (1) 

where 7 means the external electromagnetic field, t]b is the interpolating current of the 
corresponding baryon. It follows from this expression that we need the explicit expression 
of the interpolating currents to calculate the correlation function. It is well known that 
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there is a continuum number of interpolating currents for the baryons. The general form 
of the interpolating currents for the octet baryons can be represented as [IE, 
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where a, b, c are the color indices, A\ — 1, A\ — £75, A\ — 75, A\ — t and t is an arbitrary 
parameter. Ioffe current corresponds to the choice t = — 1. 

Let us firstly discuss the hadronic representation of the correlator function. By inserting 
a complete set of states between the currents in Eq. (|l|) with quantum numbers of the 
corresponding baryon B, one obtains the hadronic representation of the correlator 
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where P2 = Pi + Q, and q is the photon momentum. The second term in Eq. (|3]) takes into 
account higher states and continuum contributions and h forms a complete set of baryons 
having the same quantum number as the ground state baryon B (here B represents any 
one of p, n, S, A and S). The coupling strength of the interpolating currents with the 
baryon states parametrized with the overlap amplitudes X B is defined as 



(0\Vb\B(p)) = X B u B {p) 
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It follows from Eq. ([|) that in order to obtain an expression for the phenomenological part 
of the correlator, explicit form of the matrix element (B(pi)\B(p 2 )) 1 is needed. Electro- 
magnetic vertex of the spin 1/2 baryons can be written as 
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where the form factors fi are in general functions of q 2 = {p2—pi) 2 , and is the polarization 
vector of the photon. In the problem under consideration, only the values of the form factors 
at q 2 = are needed. 

Using Eqs. (|3])-(|5]), for the phenomenological part of the correlator ([[]) we obtain 
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where • • • represents contributions from the higher states and the continuum. As is obvious 
from Eq. (|6|), the correlator function contains some number of different structures. Among 
all possible structures, we choose the one ^$1^4 that contains the magnetic form factor 
/1 + fz- The value of j\ + f 2 at q 2 = gives the magnetic moment of the baryon in units 
of its natural magneton, i.e., eh/2rriBC. 

Isolating the structure $ 4 from the phenomenological part of the correlator, which 
describes the magnetic form factor, we get 

1 + , (7) 



Pi ~ m B P2- m B 

where /i=(f 1 + f 2 ) | ?2=0 . 

Using the explicit form of the interpolating currents given in Eq. (Q) and after some 
calculations, for the QCD part of the correlator functions of the A and S + baryons, we get 
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where C is the charge conjugation operator, subscripts a, b, c, d, e, / are the color indices and 
S q is the full propagator of the light quark involving both perturbative and nonperturbative 
parts. The expression for the theoretical parts of the correlator functions of the £ - , S°, S° 



and 5 baryons can be obtained from Eq. 
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Our calculations show that contributions from quadratic terms in strange quark mass m s 
are negligibly small compared to that of linear terms in m s , which is about 6%. Because 
of this reason theoretical part of the correlator function is calculated up to linear order in 
m s . For the full light quark propagator we have used the following expression 

S q (x) = (0\T{q(x)q(0)}\0) 

i $ m q {qq) f im q \ x 2 2/ _ . / im q 
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where A is the energy cut off separating perturbative and nonperturbative regimes. 

Few words about this expression of the light quark propagator are in order. The com- 
plete light cone expansion of the light quark propagator in external field has been carried 
in 1171 ) and it has been shown that it gets contributions from nonlocal quark operators 



qGq, qGGq, qqqq, where is the gluon field strength tensor. In the present work we 
consider the operators with only one-gluon field and neglect components with two-gluon 
and four-quark field. Formally, neglecting two-gluon and four-quark field terms can be 
justified on the basis of an expansion in conformal spin |0| . 

Perturbative part (i.e., photon interacting with quarks perturbatively) of the correlator 
function can be obtained by making the following substitution in one of the propagators in 
Eq. © 

ah 
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where the Fock-Schwinger gauge x^A^x) = is used and S^ ree = i /fc/(27r 2 x 4 ) is the 
free propagator, the remaining two propagators are the full propagators of the quarks (see 
Eq. (|9D) and J-^ v is the electromagnetic field strength tensor. The explicit expression of 
the nonperturbative contribution can be obtained from Eq. (|j) by replacing one of the 
propagators 

S q % - , (13) 

where Aj = jl, 75, 7^, 2757a, o a pj\p2^ is the full set of Dirac matrices and sum over Aj 
is implied, and the other two propagators are the full ones. 

It follows from Eqs. (|8]) and (|9|) that in order to calculate the nonperturbative contri- 
butions to the theoretical part of the correlator functions, one needs to know the matrix 
elements (7 |<?A,-g| 0) of the nonlocal operators between photon and vacuum states. Up to 
twist -4, non-zero matrix elements of the nonlocal operators are determined in terms of the 



photon wave functions as follows [[19 
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where x is the magnetic susceptibility of the quark condensate, e q is the quark charge, 
4>{u) and if)(u) are the leading twist-2 photon wave functions, while g\{u) and ^(w) are the 
twist -4 functions, respectively. Note that twist-3 photon wave functions are neglected in 
further calculations since their contribution changes the results about 5%. 

Using Eqs. (|ll]) and (0), we can calculate theoretical part of the correlator functions 
for A and S + baryons from Eqs. @ and @. The sum rule is obtained by equating the the- 
oretical and phenomenological parts of the corresponding correlator function. Performing 
double Borel transformations on the variables p\ = p 2 and p\ = (p + q) 2 in order to suppress 
the continuum and higher state contributions (for a discussion concerning this point, see 
21 — 1|23||)) the following sum rules for the A and E + baryons are obtained: 
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are used to subtract the continuum and higher state contributions, x = s /M 2 and s is 
the continuum threshold. Moreover 
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but, since we are dealing with just a single baryon, the Borel parameters Mf and Mf should 
be set to equal each other, from which it follows that Uq — 1/2. 

As has already been noted, the sum rules for the E°, £~, H°, S~, p, n baryons can be 
obtained by making the following replacements: 
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Note that, after we make these replacements we set masses of u and d quarks to zero and 
assume SU{2) flavor symmetry, which implies that (mm) = (dd). 

Obviously, it follows from Eqs. (|TJ) and flTBD that one needs to know the residues A# 
of the octet baryons in order to determine their magnetic moments. These residues are 
determined from mass sum rules of the corresponding baryons: 
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3 Numerical Analysis 

In this section we present the numerical analysis of the sum rules for the magnetic moments 
of the octet baryons. The main input parameters in the LCQCD are the photon wave 
functions which have been introduced in the previous section. In [18] and it has been 



shown that the leading photon wave functions do not deviate much from their asymptotic 
form. Therefore in numerical calculations we shall use the following forms of the photon 
wave functions iI8L ED 



4>{u) = 6u(l — u) , ip(u) = 1 , 

I 1 
9i(u) = --(1 - u)(3 - u) , g 2 (u) = --(1 - u) 2 . 

The values of the other input parameters that are used in the numerical analysis are: / = 
0.028 GeV 2 , x = -3.3 GeV~ 2 @ (in this quantity is estimated to be % = -4.4 GeV~ 2 ^ 



{qq)(l GeV) = -(0.243) 3 GeV 3 , m 2 = (0.8±0.2) GeV 2 g§, m s (l GeV) = (150±50) MeV, 
(ss)(l GeV) = 0.8(gg)(l GeV) and A = 0.5 GeV. In the problem under consideration we 
have three auxiliary parameters, namely parameter t in the hadron interpolating currents, 
Borel mass square M 2 and the continuum threshold so- No doubt, it is expected that any 
physically measurable quantity must be independent of these auxiliary parameters. So, our 
problem is to find the appropriate regions, for which magnetic moments of octet baryons 
are independent of the above-mentioned parameters. 

For this aim we consider the following three-step procedure. In the first step we attempt 
to find working region of M 2 , where magnetic moments of octet baryons are independent of 
the Borel parameter at fixed values of Sq and t. In Figs. (l)-(8) we present the dependence 
of the magnetic moments of the corresponding baryons on M 2 . From these figures we 
see that the magnetic moments of octet baryons seem to be almost independent of M 2 
for different choices of t and sq. However, the working region for the Borel mass square 
parameter for the members of octet baryons are different, i.e., 0.9 GeV 2 < M 2 < 1.2 GeV 2 
for p, n; 1.3 GeV 2 < M 2 < 1.6 GeV 2 for E~, S°, £+ and A; 1.7 GeV 2 < M 2 < 2.1 GeV 2 
for S° and H~. 

Before determining the magnetic moment of baryons, the next problem to be considered 
is to find the working region of t, i.e. the region where the sum rules make sense. For this 
purpose we consider the mass sum rules (Eqs. ([17|-p0D). Two criteria should be satisfied 



for mass sum rules. First of all, each mass sum rule must separately be positive. The 
second criteria is that, the predicted mass of the baryons which is obtained by taking 
the logarithmic derivative of the corresponding sum rule with respect to M~ 2 should be 
stable with respect to the variation of the parameter t. As a result of analysis of the mass 
sum rules we found that, the region —0.5 < t < 0.5 is unphysical (see Fig. (9)). Our 
analysis shows that the working region of t where the second criteria is satisfied is given by 
—0.6 < cos# < 0.3, where 9 is defined through the relation 9 = tan^i. This region of 9 
corresponds to t < —1.4 and t > 3.3. We observe that the Ioffe current, which corresponds 
to the choice t = — 1 does not lie in the working region. 
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Our last attempt is to determine a region for the parameter t where the magnetic 
moments of the octet baryons are independent of its variation. Here we stress again that 
in the first step the working region for the Borel mass square M 2 has been determined for 
which [Ib is practically independent of its variation, as well as insensitive to the continuum 
threshold. We present in Fig (10) the dependence of the magnetic moments of the octet 
baryons (in units of the nucleon magneton fi^) on cos 6, at fixed values of M 2 (remember 
that the fixed values of M 2 are different for each one of the octet baryons) and at fixed value 
of So (again different fixed values of Sq are used for each baryon). From this figure we observe 
that hb is practically independent of the parameter t in the region —0.6 < cos 6 < 0.3 
and insensitive to the variation of the continuum threshold sq. Our final results on the 
magnetic moments of the octet baryons are presented in Table (1) . For completeness, in 
this table we also present the predictions of other approaches on the magnetic moments 
of octet baryons and the experimental results as well as our results for the value of the 
magnetic susceptibility x = — 4.4 GeV 2 . It follows from our analysis that light cone QCD 
sum rule prediction on magnetic moments of octet baryons are close to their experimental 
values for x — —3.3 GeV 2 , but considerable departure is observed for x — —4.4 GeV 2 . 
Obviously, it follows from Table 1 that similar situation takes place for the traditional sum 
rule predictions on the magnetic moments of octet baryons. The uncertainties in the results 
we have presented could be attributed to the variation of so, Borel parameter M 2 , twist -3 
photon wave functions which are neglected and errors in the input parameters such as (qq), 
t«q and m s . All errors are added quadratically to the final predictions. It should be stressed 
here that the predictions on magnetic moments at t = — 1, which corresponds to the Ioffe 
current, are not reliable. 

From our results we see that the values of ~ 2/ia ~ /i n ; Hs° ~ and ~ /is- . 
It should be added here that, in exact SU(3) symmetry, the relations between the magnetic 



moments £1X6 cLS follows |3 



/^E+ — Hp > 
2/iA = fi n , 

fj, s - = /Is- , 

/j, s o = fj, n . (21) 

Violation to the relations ( pip comes from SU(3) breaking terms (the mass of the s quark, 
s quark condensate, etc.). So, our predictions give us a clue about sizable SU(3) symmetry 
breaking effects. 

In summary, we have calculated magnetic moments of the octet baryons in light cone 
QCD sum rules, using the general form of the baryonic interpolating currents. We have 
obtained that light cone QCD sum rule predictions on magnetic moments of the octet 
baryons are in good agreement with the experimental data at x = — 3.3 GeV 2 . 
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Note added: After we have completed the present work, we have received a recent work 
where the mass sum rules for general form of the interpolating currents for the baryons 



are calculated. Our results on mass sum rules coincide with with the ones predicted in this 
work. 
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NQM 


SQM 


QCDSR 


QCDSA 




SKRM 


LCQSR 


EXP 


X = -3.3 


X = -4-4 


X = -3-3 


X = -4-4 


pp 


9 87 


2 75 


z. / z 


O.00 


2 ^4 


2 793 


2 36 


z. ( ± U.o 


O.O ± U.O 


2 79 




—1 91 


— 1 84 


— 1 65 


—2 06 

Zj . V V .' 


— 1 69 


— 1 913 


— 1 87 


— 18 + 35 

-L.O _l_ V J . (J '.J 


—2.3 ± 0.4 


— 1 91 


1 z_i 


2.62 


2.65 


2.52 


3.30 


2.48 


2.458 


2.46 


2.2 ±0.4 


2.9 ±0.4 


2.46 ±0.01 




-1.20 


-1.02 


-1.13 


-1.38 


-0.90 


-1.16 


-1.16 


-0.8 ±0.2 


-1.1 ±0.3 


-1.16 ±0.03 




-0.63 


-1.44 


-0.89 


-0.98 


-1.49 


-1.25 


-1.25 


-1.3 ±0.3 


-1.3 ±0.4 


-1.25 ±0.01 


Hr- 


-1.44 


-0.52 


-1.18 


-1.27 


-0.63 


-0.6531 


-0.65 


-0.7 ±0.2 


-1.0 ±0.2 


-0.65 




-0.63 


-0.67 


-0.50 


-0.80 


-0.69 


-0.613 


-0.60 


-0.7 ±0.2 


-0.9 ±0.2 


-0.61 



Table 1: Predictions of various approaches for the octet baryon magnetic moments: naive quark model (NQM, see ref. in [31]); 
static quark model (SQM) [32]; QCD sum rules (QCDSR) [9]; QCD string approach (QCDSA) [28]; chiral perturbation theory 
(xPT) [29]; skyrme model (SKRM) [30]; present work (LCQSR). For completeness we present the experimental values of the 
octet baryons. All the values in the table are given in units of nuclear magneton fi^. 
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Figure captions 



Fig. (1) The dependence of the magnetic moment (in units of nuclear magneton /in) of p 
on the Borel parameter M 2 , at the continuum threshold s = 2.0 GeV 2 . 

Fig. (2) The same as in Fig. (1), but for n. 

Fig. (3) The same as in Fig. (1), but for A, at the continuum threshold s = 2.5 GeV 2 
and m s = 0.15 GeV. 

Fig. (4) The same as in Fig. (3), but for S°, at the continuum threshold s = 3.0 GeV 2 
and m s = 0.15 GeV. 

Fig. (5) The same as in Fig. (4), but for £ _ . 
Fig. (6) The same as in Fig. (4), but for S + . 

Fig. (7) The same as in Fig. (3), but for S°, at the continuum threshold s — 3.2 GeV 2 
and m s = 0.15 GeV. 

Fig. (8) The same as in Fig. (7), but for S _ . 

Fig. (9) The dependence of \ 2 B on t, for all members of the octet baryons, at the contin- 
uum threshold s = 2.0 GeV 2 (for p and n),s = 2.5 GeV 2 (for A), s — 3.0 GeV 2 (for £) 
and s = 3.2 GeV 2 (for S). 

Fig. (10) The dependence of the magnetic moments of the octet baryons on cos# at the 
continuum threshold So = 2.0 GeV 2 (for p and n), s = 2.5 GeV 2 (for A), s = 3.0 GeV 2 
(for S) and s = 3.2 GeV 2 (for S), and at M 2 = 1 Ge\/ 2 (for p and n), M 2 = 1.5 GeV 2 
(for E and A) and M 2 = 1.9 Ge\/ 2 (for S). 
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